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Abstract The study of the metal enrichment of the intra-cluster and inter-galactic media 
(ICM and IGM) represents a direct means to reconstruct the past history of star formation, 
the role of feedback processes and the gas-dynamical processes which determine the evolu- 
tion of the cosmic baryons. In this paper we review the approaches that have been followed 
so far to model the enrichment of the ICM in a cosmological context. While our presenta- 
tion will be focused on the role played by hydrodynamical simulations, we will also discuss 
other approaches based on semi-analytical models of galaxy formation, also critically dis- 
cussing pros and cons of the different methods. We will first review the concept of the model 
of chemical evolution to be implemented in any chemo-dynamical description. We will em- 
phasise how the predictions of this model critically depend on the choice of the stellar initial 
mass function, on the stellar life-times and on the stellar yields. We will then overview the 
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comparisons presented so far between X-ray observations of the ICM enrichment and model 
predictions. We will show how the most recent chemo-dynamical models are able to cap- 
ture the basic features of the observed metal content of the ICM and its evolution. We will 
conclude by highlighting the open questions in this study and the direction of improvements 
for cosmological chemo-dynamical models of the next generation. 

Keywords Cosmology: numerical simulations - galaxies: clusters - hydrodynamics - 
X-ray: galaxies 



1 Introduction 

Clusters of galaxies are the ideal cosmological signposts to trace the past history of the 
inter-galactic medium (IGM), thanks to the high density an d temperature reached by the 
cosmic baryons trap ped in their gravitational potential wells feosati et alJl200i |wi|[2005 ; 
iDiaferio et ail [2008 - Chapter 2, this volume). Observations in the X-ray band with the 
Chandra and XMM-Newton satellites are providing inval uable information on the thermo- 
dynamical properties of the intra-cluster medium (ICM) ( iKaastra et al.ll2008l - Chapter 9, 
this volume). These observations highlight that non-gravitational sources of energy, such as 
energy feedback from supernovae (SNe) and Active Galactic Nuclei (AGN) have played an 
important role in determining the ICM physical properties. Spatially resolved X-ray spec- 
troscopy permits to measure the equivalent width of emission lines associated to transi- 
tions of heavily ionised elements and, therefore, to trace the pattern of chemical enrichment 
fe.g.. lMushotzkvll2004l for a review). In turn, this information is inex tricably lin ked to the 
history of formation and evolution of the galaxy population (e.g.. lRenzinilll997l and ref- 
erences therein) , as in ferred from observations in the optical/near-IR band. For instance, 
be Grandi et al.1 d2004l ) have first shown that cool core clusters are characterised by a sig- 
nificant central enhancement of the iron abundance, which closely correlates with the mag- 
nitude of the Brightest Cluster Galaxies (BCGs). This demonstrates that a fully coherent 
description of the evolution of cosmic baryons in the condensed stellar phase and in the dif- 
fuse hot phase requires properly accounting for the mechanisms of production and release 
of both energy and metals. 

The study of how these processes take place during the hierarchical build up of cosmic 
structures has been tackled with different approaches. Semi-analytical models (SAMs) of 
galaxy formation provide a flexible tool to explore the space of parameters which describe 
a number of dynamical and astrophysical processes. In their most recent formulation, such 
models are coupled to dark matter (DM) cosmological simulations, to trace the merging 
history of the halos where galaxy formation takes place, and include a tre atment of metal 
production from type-la and type-II supernovae (SNIa and SNII, hereafter: iDe Lucia et al.l 
l2004l : lNagashima et alj2005h . so as to properly address the study of the chemical enrichment 
of the ICM. The main limitation of this method is that it does not include the gas dynamical 
processes which causes metals, once produced, to be transported in the ICM. As a conse- 
quence, they provide a description of the global metallicity of clusters and its evolution, but 
not of the details of its spatial distribution. 

In order to overcome this limitation, ICoral J2006l) applied an alternative approach, in 
which non-radiative hydrodynamical simulations of galaxy clusters are used to trace at the 
same time the formation history of DM halos and the dynamics of the gas. In this approach, 
metals are produced by SAM galaxies and then suitabl y assigned to gas parti cles, thereby 
providing a chemo-dynamical description of the ICM. Domainko et al.l j2006l) used hydro- 
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dynamical simulations, which include prescriptions for gas cooling, star formation and feed- 
back, similar to those applied in S AMs, to address the sp ecific role played by ram-pressure 
stripping to distribute metals, while lKapferer et all d2007bh used the same approach to study 
the different roles played by galactic winds and by ram-pressure stripping. While these ap- 
proaches offer advantages with respect to standard SAMs, they still do not provide a fully 
self-consistent picture, in which chemical enrichment is the outcome of the process of star 
formation, associated to the cooling of the gas infalling in high density regions. In this sense, 
a fully self-consistent approach requires that the simulations include the processes of gas 
cooling, star formation and evolution, along with the corresponding feedback in energy and 
metals. 

A number of authors have presented hydrodynamical simulations for the formation of 
cosmic stru ctures, which includ e treatments of the chemical evolution at different levels of 
complexity. iRaiteri et al.l ll 19961) presented SPH simulations of the Galaxy, forming in an 
isolated halo, by following iron and oxygen production fro m SN II and SN la, also account- 
ing for the effect of stellar lifetimes. iMosconi et al.l d200lh performed a detailed analysis of 
chemo-dynamical SPH simulations, aimed at studying both numerical stabi lity of the results 
and the enrichment properties of galactic objects in a cosmological context- Ilia et al. 
discussed a statistical approach to follow metal production in SPH simulations, which have a 
large number of sta r particles, showing applica tions to simulations of a disc-like galaxy and 
of a galaxy cluster. iKawata & Gibsonl j2003bh carried out cosmological chemo-dynamical 
simulations of elliptical galaxies, with an SPH code, by including the contribution from 
SNIa and SNII, also accounting for stellar lifetimes. IValdamirul d2003h performed an ex- 
tended set of cluster simulations and showed that profiles of the i ron abundanc e are s teeper 
than the observed ones. A similar analysis has been presented by Ro meo et al.l j2006h . who 
also considered the effe ct of varying the IMF an d the feedback efficiency on the enrich- 
ment pattern of the ICM. Scannapieco et al. (2005) presented an implementation of a model 
of chemical enrichment in the GADGET-2 code, coupled to a self-consistent model for star 
formation and feedback. In their model, which was applied to study the enrichment of galax- 
ies, they included the contri bution from SNIa and SN II, assuming that all long-lived stars 
die after a fixed delay time. iTornatore et al. U2007al . l2004) presented results from an imple- 
mentation of a detailed model of chemical evolution model in the GADGET-2 code dSpringell 
l2005h . including metallicity-dependent yields and the contribution from intermediate and 
low mass stars (ILMS hereafter). The major advantage of this approach is that the metal 
production is self-consistently predicted from the rate of gas cooling treated by the hydro- 
dynamical simulations, without resorting to any external model. However, at present the 
physical scales involved by the processes of star formation and SN explosions are far from 
being resolved in simulations which sample cosmological scales. For this reason, such sim- 
ulations also need to resort to external sub-resolution models, which provide an effective 
description of a number of relevant astrophysical processes. 

The aim of this paper is to provide a review of the results obtained so far in the study 
of the chemical enrichment of the ICM in a cosmological context. Although we will con- 
centrate the discussion on the results obtained from full hydrodynamical simulations, we 
will also present re sults based on SAMs. As such, this paper complements the reviews by 
IBorgani et ai]|2008l - Chapter 13, this volume, which reviews the current status in the study 
of the thermod ynamical properties of the ICM with cosmological hydrodynamical simu- 
lations, and by ISchindler & Diaferidl2008l - Chapter 17, this volume, which will focus on 
the study of the role played by different physical processes in determining the ICM enrich- 
ment pattern. Also, this paper will not present a detailed description of the techniques for 
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simulations of galaxy clusters, which is reviewed by iDolag et alJEooi - Chapter 12, this 
volume. 

In Sect.[2]we review the concept of model of chemical evolution and highlight the main 
quantities which are required to fully specify this model. In Sect. [3] we review the results 
on the global metal content of the ICM, while Sect. |4] concentrates on the study of the 
metallicity profiles and Sect. [5] on the study of the ICM evolution. Sect. [6] discusses the 
properties of the galaxy population. Finally, Sect.|7]provides a critical summary of the results 
presented, by highlighting the open problems and lines of developments to be followed by 
simulations of the next generation. 



2 What is a model of chemical evolution? 

In this section we provide a basic description of the key ingredients required b y a model 
of che mical evolution. For a more detailed description we refer to the book by iMatteuccl 
d2003h . 

The process of star formation in cosmological hydrodynamical simulations is described 
through the conversion of a gas element into a star particle; as a consequence, in simulations 
of large volumes the star particles have a mass far larger than that of a single star, with values 
of the order of 10 6 — 10 8 M ffi , de pending on the resolution and on the mass of the simulated 
structure (e.g. jKatz et al ] | 19961) . The consequence of this coarse-grained representation of 
star formation is that each star particle must be treated as a simple stellar population (SSP), 
i.e. as an ensemble of coeval stars having the same initial metallicity. Every star particle 
carries all the physical information (e.g. birth time t^, initial metallicity and mass) that is 
needed to calculate the evolution of the stellar populations that they represent, once the 
lifetime function (see Sect. |2.2| |, the IMF (see Sect. |2.4| | and the yields (see Sect. |2.3I > for 
SNe and ILMS have been specified. Therefore, we can compute for every star particle at any 
given time t > fb how many stars are dying as SN II and SN la, and how many stars undergo 
the AGB phase, according to the equations of chemical evolution, that we will present here 
below. 

The generation of star particles is pertinent to hydrodynamical simulations, which in- 
clude the description of the star formation. As for semi-analytic models (SAMs) of galaxy 
formation, they generally identify a galaxy with a suitably chosen DM particle in the col- 
lisionless simulation which is used to reconstruct the merger tree. While this DM particle 
defines the position of a galaxy, the corresponding stellar content is described as a superpo- 
sition of SSPs, each generated i n correspondence of t he time step with which the process of 
galax y formation is studied (e.g jDiaferio et al.l200ll : IBe Lucia et al.l2004l :lNagashima et alJ 
120051 " for applications of SAMs to the cluster galaxies). In this sense, the description of the 
chemical evolution model, as we provide here, can be used for both hydrodynamical simu- 
lations and for SAMs. 

It is generally assumed that the stars having mass above 8 M Q at the end of the hydro- 
static core burning undergo an electron capture process, leading to a core collapse. A large 
amount of energy can be transferred to the outer layers during this phase due to a substantial 
production of neutrinos that easily escape from the central core. Although theoretical work 
has not yet been able to reproduce a sufficient energy deposition, it is currently supposed 
that this process leads to an explosive ejection of the outer layers, giving rise to a SN II. We 
remind the reader that 8 M a is a commonly adopted fiducia l value, although the limiting 
mass for the onset of explosive evolution is still debated (e.g.. |Portinari et al. [l998). 
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A different ejection channel is provided by the SN la that are believed to arise from ther- 
monuclear explosions of white dwarfs (WD hereafter) in binary stellar system s as a con- 
sequence of the matter accretion from the companion (e.g.. lNomoto et ai] |2000). However, 
there are still a number of uncertainties about the nature of both the WD and the companion 
and about the mass reach ed at the onset of the explosion (e.g.. iMatteucci & Recc h1 l200lL 
lYungelson & Lividfc oOO ) . Furthermore, observational evidences for multiple populations of 
Type la progenitors have been recently found both from direct detection of Ty pe la events 

200 6|) or from in 



(e.g.. lMannucci et al.ll2005l . lScannapieco & Bildstenll2005l , ISullivan et al 



feren ces on the chemical enrichment patterns in galaxy clusters (e.g., I Gal- Yam & Maozl 
120041) . Finally, a third way to eject heavy elements in the interstellar medium is the mass 
loss of ILMS by stellar winds. 

In summary, the main ingredients that define a model of chemical evolution are the 
following: (a) the SNe explosion rates, (b) the adopted lifetime function, (c) the adopted 
yields and (d) the IMF which fixes the number of stars of a given mass. We describe each of 
these ingredients in the following. 



2.1 The equations of chemical evolution 
2.7.7 Type la supernovae 

We provide here below a short description of how Type la SN are included in a model of 
chemical evolution. For a comprehe nsive review of analytical formulations we refer to the 
paper by Greggio (2005). Following IGreggio & Renzinil dl983l ). we assume here that SNIa 
arise from stars belonging to binary systems, having mass in the range 0.8-8 M . Accord- 
ingly, in the single-degenerate WD scenario dNomoto et alj|2000h . the rate of explosions of 
SN la can be written as 

.sup Mm 

7?SNIa(0= A J 4>('«b) / /(/i)V / '( f -fm 2 )dM dm B. (1) 
M B .M Mm 

In the above equation, (j)(m) is the IMF, is the total mass of the binary system, m% is 
the mass of the secondary companion, x m is the mass-dependent life-time and \j/(t) is the 
star formation rate. The variable ji = mi/ma is distributed according to the function /(/x)> 
while A is th e fraction of stars in binary systems of that particular type to be progenitors 



in binary 

of SNIa (see Matteucci & Recchi 2001 for more details). For instance, in the model by 
IGreggio & RenziniT i 19831) fj, varies in the range between fi m and fi M = 0.5, with ji„ 



max [mzft) /ms, (ff? B — 0.5M B m)/'« B ], where triiii) is the mass of the companion which dies 
at the time t, according to the chosen life-time function. Furthermore, let Msm an d M B m be 
the smallest and largest value allowed for the progenitor binary mass m B . Then, the integral 
over ;«b runs in the range between M B m f and M Bsup , which represent the minimum and 
the maximum value of the total mass of the binary system that is allowed to explode at 
the time t. These values in general are functions of M Bm , M B m, and mi(f), which in turn 
depends on the star formation history f(f). The exact func tional d ependence is defined by 
the SNIa progenitor model. For instance, in the model by Greg gio & Renzinil (l983) it ls 
^B.inf = max[2m2(/),M Bm ] and M B sup = 0.5M BM +m2(f). Under the assumption of a short 
duration burst of star formation, the function \j/(t) can be approximated with a Dirac 5- 
function. This case applies to hydrodynamical simulations, which include star formation, 
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where the creation of a SSP is described by an impulsive star formation event, while more 
complex descriptions should take into account the continuous star formation history yf(t). 
Under the above assumption for \j/(t) and using the functional form of f(n) derived from 
statis t ical studies of the stellar population in the solar neighbourhood ( ITutukov & Iungelsonl 
ll98dl : lMatteucci & RecclrfeOOll) . we find 



flSNIaW 



dm 2 {t) 



dt 



m 2 =T~ [ (t) 



24 mj A 



Mbm 1 

(j)(m B )—rdm B . 



(2) 



Since the current understanding of the process of star formation does not allow to precisely 
determine the value of A, its choice can be fixed from the requirem ent of reproducing a 
specifi c observation, once the form of the IMF is fixed. For instance, Matt eucci & Gibsonl 
(1995) found that A = 0.1 was required to reproduce the observed iron enrichment. 

2.7.2 Supernova Type II and low and intermediate mass stars 

Computing the rates of SN II and ILMS is conceptually simpler, since they are driven by the 
lifetime function %{m) convolved with the star formation history \j/(t) and multiplied by the 
IMF <t>(m = T (f)). Again, since \j/(t) is a delta-function for the SSP used in simulations, 
the SN II and ILMS rates read 



R 



SNII\UMS 



(t) = (t>(m{t))x 



dm(t) 
dt 



(3) 



where m(t) is the mass of the star that dies at time t. We note that the above expression must 
be multiplied by a factor of (I — A) for AGB rates if the interested mass m(t) falls in the 
same range of masses which is relevant for the secondary stars of SNIa binary systems. 

In order to compute the metal release by stars (binary systems in case of SN la) of a 
given mass we need to take into account the yields pzi(m,Z), which provide the mass of the 
element i produced by a star of mass m and initial metallicity Z. Then, the equation which 
describes the evolution of the mass p, (f) for the element holding for a generic form of the 
star formation history \jf(t), reads: 

p,(r) = -w{t)Zi{t) 



+ A J <j){m) 



0.5 

/ f(LiW(t-T m2 )pz,(m,Z)dp: 



dm 



Mbm 

+ (1-A) / \jf(t - x(m))pz i (m,Z)<p(m)dm 



(4) 



A^Bir 



+ / \jf(t - l{m))p Zi {m 1 Z)<p{m)dm 



+ J \jf(t— T(m))pz i (m,Z)(p(m)dm. 

In the above equation, Ml and M\j are the minimum and maximum mass of a star, re- 
spectively. Commonly adopted choices for these limiting masses are Ml ^ 0.1 M a and 
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M\j ~ 100 M Q . The term in the first line of Eq.|4]accounts for the metallicity sink due to the 
locking of metals in the newborn stars. The term in the second line accounts for metal ejec- 
tion contributed by SN la. The terms in the third and fourth lines describe the enrichment by 
mass-loss from intermediate and low mass stars, while the last line accounts for ejecta by 
SNH 



2.2 The lifetime function 



Different choices for the mas s-dependence of the l i fe-tim e function have been proposed 
in the literature. For instance, Padovani & M atteuccil d 19931) (PM93 hereafter) proposed the 



expression 



10 [(l-34-Vl-79-°-22(7.76-log(m)))/0.11]-9 for m < g 6 M 



x(m) 



(5) 



1.2m 



-1.85 



+ 0.003 otherwise. 



An alternative expression has been propos e d by Maeder & Mevnetl dl989l) (MM89 here- 
after), and extrapolated by IChiappini et all d 19971) 



to very high (> 60 M ) and very low 



(< 1.3 M Q ) masses: 



T(m) : 



1Q-O.65451ogm+1 
jQ-3.71ogm+1.351 
lQ-2.511ogm+0.77 
]Q-1.781ogm+0.17 
jQ-0.S61ogm-0.94 

1.2 x m~ 185 +0.003 otherwise 



m< 1.3 M Q 
1.3<m<3M 
3 < m < 7 Mro 
7 <m< 15 M Q 
15 < m < 53 M 



(6) 



The main difference between these two functions concerns the life-time of low mass 
stars (< 8M ra ). The MM89 life-time function delays the explosion of stars with mass > 
1 M a , while it anticipates the explosion of stars below 1 M a with respect to the PM93 
life-time function. Only for masses below 1 M Q does the PM93 function predict much 
more long-living stars. This implies that different life-t imes will produce diff erent evolution 
of both absolute and relative abundances (we refer to iRomano et alJ d2005l) for a detailed 
description of the effect of the lifetime function in models of chemical evolution). 

We point out that the above lifetime functions are independent of metallicit y, whereas in 
princ iple this depende nce can be included in a model of chemical evolution (e.g. JPortinari et al] 
Il9980 . For instance, Rai teri et alJ dl99q) u sed the metallicity -dependent lifetimes as ob- 
tained from the Padova evolutionary tracks ( IBertelli et alJll994l ). We show in Fig.[T]a com- 
parison of different lifetime functions. 



2.3 Stellar yields 

The stellar yields specify the quantity pZj{m,Z), which appears in Eq. |4]and, therefore, 
the amount of different metal species which are released during the evolution of a SSP 
A nu mber of diffe r ent se ts of yields have been proposed in the l iterature , such as those 
by iRenzini & VbiH dl98ll) . Ivan den Hoek & Groenewegenl d 19971) . iMarigol d200ll) for the 
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Fig. 1 Left panel: the mass dependence of different lifetime functions: Padovani & Matteucci (1993, solid 
line), Maeder & Meynet (1998, dotted line), Raiteri et al. (1996, dashed lines; the thin line refers to Z = 
O.O2Z0, the thick line refers to Z = 2Zq), Portinari et al. (1998, dot-dashed lines; the thin line refers to 
Z = O.O2Z0, the thick line refers to Z = 2Z R ). Right panel: the ratio between the same lifetimes and that by 
Padovani & Matteucci (1993). 



LIMS and those bv lNomoto etail d 19971) . Ilwamoto et alj dl999l) . iThielemann et alj [20031) 
for SN la. As fo r SNII, there are many prop ose d sets of metallicity-d epe ndent yields; among 
others, those bv lWooslev & Weaver! dl995l) . bv lPortinari et all dl998l) . bv lChieffi & Limongil 
d2004l) . which are based on different assumptions of the underlying model of stellar structure 
and evolution. 

As an example of the differences among different sets of yields, we show in the left panel 
of Fig. [2] the ratios betwe en the abundances o f diffe rent elemen ts produced by the S NII o f 
a SSP, as expected from IWooslev & Weaver! d 19951) and from IChieffi & Limongil d2004l) . 
Different curves and symbols here correspond to different values of the initial metallicity of 
the SSP. Quite apparently, the two sets of yields provide significantly different metal masses, 
by an amount which can sensitively change with initial metallicity. This illustrates how a 
substantial uncertainty exists nowadays about the amount of metals produced by different 
stellar populations. There is no doubt that these differences between different sets of yields 
represent one of the main uncertainties in any modelling of the chemical evolution of the 
ICM. 

For this reason, we recommend that papers related to the cosmological modelling of 
the chemical evolution, either by SAMs or by hydrodynamical simulations, should be very 
careful in specifying for which set of yields the computations have been performed as well 
as the details of assumed lifetime and models for SNII and SNIa. In the absence of this, 
it becomes quite hard to judge the reliability of any detailed comparison with observational 
data or with the predictions of other models. 



2.4 The initial mass function 

The initial mass function (IMF) is one of the most important quantities in a model of chemi- 
cal evolution. It directly determines the relative ratio between SN II and SN la and, therefore, 
the relative abundance of a -elements and Fe-peak elements. The shape of the IMF also de- 
termines how many long-living stars will form with respect to massive short-living stars. In 
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turn, this ratio affects the amount of energy released by SNe and the present luminosity of 
galaxies, which is dominated by low mass stars, and the (metal) mass-locking in the stellar 
phase. 

As of today, no general consensus has been reached on whether the IMF at a given time 
is universal or strongly dependent on the environment, or whether it is time-dependent, 
i.e. whether local variations of the values of temperature, pressure and metallicity in star- 
forming regions affect the mass distribution of stars. 

The IMF (j> (m) is defined as the number of stars of a given mass per unit logarithmic 
mass interval. A widely used form is 

<j)(m) = dN/dlogm °c m~ x{ - m) . (7) 

If the exponent x in the above expression does not depend on the mass m, the IMF is then 
descri bed by a single power-law. The most fa mous and widely u s ed sin gle power-law IMF 
is the lSalpeteJ dl955b one that has x = 1.35. lArimoto & Yoshiil l l 19871 . AY hereafter) pro- 
posed an IMF with x = 0.95, which predicts a relatively larger number of massive stars. In 
general, IMFs providing a large number of massive stars are usually called top-heavy. More 
recently, different expressions of the IMF have been proposed in order to model a flatten- 
ing in the low-mass regime that is currently favoured by a number of observations. IKroupal 
fcooih introduced a multi-slope IMF, which is defined as 



(j)(m) <x < 



m~ 13 ?n>0.5M e 

m-° 3 0.08<m<0.5M e (8) 
m 01 m< 0.08 M m 



IChabried J2003h proposed another expression for the IMF, which is quite similar to that one 
proposed by Kroupa 

!m~ 13 m > 1 Mq 

Theoretical arguments (e.g. lLarson|[l998l) suggest that the present-day characteristic mass 
scale, where the IMF changes its slope, ~ 1 M a should have been larger in the past, so that 
the IMF at higher redshift was top-heavier than at present. 

While the shape of the IMF is determined by the local conditions of the inter-stellar 
medium, direct hydrodynamical simulations of star formation in molecular clouds are only 
now approaching the required resolution and sophistication le vel to make credible predic- 
tions on the IMF (e.g jBonnell et al.ll2006l . |Padoan et al.ll200l see lMcKee & Ostrikei|[2007l 
for a detailed discussion). 

We show in Fig.[2]the number of stars, as a function of their mass, predicted by different 
IMFs, relative to those of the Salpeter IMF. As expected, the AY IMF predicts a larger 
number of high-mass stars and, correspondingly, a smaller number of low-mass stars, the 
crossover taking place at ~ 2 M Q . As a result, we expect that the enrichment pattern of the 
AY IMF will be characterised by a higher abundance of those elements, like oxygen, which 
are mostly produced by SN II. Both the Kroupa and the Chabrier IMFs are characterised by a 
relative deficit of very low-mass stars and a mild overabundance of massive stars and ILMS. 
Correspondingly, an enhanced enrichment in both Fe-peak and a elements like oxygen is 
expected, mostly due to the lower fraction of ma ss locked in ever-liv ing stars. For reference 
we also show a different IMF, also proposed bv lKroupa etaH dl993l) , that exhibits a deficit 
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Fig. 2 Left panel: the ratio Mj v w /A^' L between the mass of species ;', produced by the SNII of a SSP, when 
using the two sets of yields by Wooslev & Weaver 1 1995) and by Chieffi & Limongi (2004) for different 
values of the initial SSP metallicity. Dif ferent symbols are for d ifferent values of the initial metallicity of 
the SSP, as reported by the labels. From Tornatore et al. (2007a). Right panel: the shape of different IMFs 
introduced in the literature. Stellar masses are expressed in units of M . 



in both very low- and high-mass stars; for this kind of IMF a lower a / Fe ratio with respect 
to the Salpeter IMF is expected. 

Since clusters of galaxies basically behave like "closed boxes", the overall level of en- 
richment and the relative abundances should directly reflect the condition of star formation. 
While a number of studies have been performed so far to infer the IMF shape fro m the en- 
richment pat tern of the ICM, no gen eral consensus has been reached. For instance, iRenzinil 
dl997l . l2004) and lPipino eTal] d2002h argued that both the global level of ICM enrichment 
and the a/Fe relative abundance can be reproduced by assuming a Salpeter IMF, as long as 
this relative abundance is sub-solar. Indeed, a different conclusion has been reached by other 
authors in the attempt o f explaining the relativ e abundances [a/Fe]^ in the ICM (e.g., 
iBaumgartner et al1l2005h - IPortinari et al.l j2004l) used a phenomenological model to argue 
that a standard Salpeter IMF cannot account for the observed a/Fe ratio i n the ICM. As we 
shall d iscuss in the following, a similar conclusion was also reached by iNagashima et alj 
J2005h . who used semi -analytical models of galaxy formation to trace the production of 



heavy elements, and by Romeo et al 



! 2005h . who used hydrodynamical simulations includ- 



ing chemical enrichment. ISaro et al.l 1 120060 analysed the galaxy population from simulations 
of galaxy clusters and concluded that a Salpeter IMF produces a colour-magnitude relation 
that, with the exception of the BCGs, is in reasonable agreement with observations. On the 
contrary, the stronger enrichment provided by a top-heavier IMF turns into too red galaxy 
colours. 

In the following sections, the above results on the cosmological modelling of the ICM 
will be reviewed in more detailed, also critically discussing the possible presence of obser- 
vational biases, which may alter the determination of relative abundances and, therefore, the 
inference of the IMF shape. 

Different variants of the chemical evolution model have been implemented by differ- 
ent authors in their simulation codes. For instance, the above described model of chem - 
ical evolution has been implemented with minimal variants bv lKawata & Gibs"or] 120 03a) 



Kobavashi ( 2004), who also included the effect of hypernova explosions, and Tornatore et alj 
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(2007a), while iTomatore et al 



(2007b) also included the effect of metal enrichment from 



low-metallicity (Pop III) stars. iRaiteri et al.ldl996t) and lValdarninil J2003h als o used a similar 
model , but neglected the contribution from low- and intermediate-mass stars. iMosconi et all 
J200ll) and Scannap ieco et al.l ( 2005) neglected delay times for SNII, assumed a fixed delay 
time for SNIa and neglected the contribution to enrichment from low- and intermediate- 
mass stars. 

Clearly, a delicate point in hydrodynamical simulations is deciding how metals are dis- 
tributed to the gas surrounding the star particles. The physical mechanisms actually respon- 
sible for enriching the inter-stellar medium (ISM; e.g., stellar winds, blast waves from SN 
explosions, etc.) take place on scales which are gener ally well below the resoluti on of cur- 
rent cosmological simulations of galaxy clusters (see [Schindler & Diaferioll2008l - Chapter 
17, this volume). For this reason, the usually adopted procedure is that of distributing met- 
als according to the same kernel which is used f or the computation of the hydrodynamical 
forces, a choice which is anyway quite arbitrary. IMosconi et al] J200ll) and lTornatore et all 
(2007a) have tested the effect of changing in different ways the weighting scheme to dis- 
tribute metals and found that final results on the metal distribution are generally rather sta- 
ble. Although this result is somewhat reassuring, it is clear that this warning on the details 
of metal distribution should always be kept in mind, at least until our simulations will have 
enough reso lution and accurate descrip tion of the physical processes determining the ISM 
enrichment dSchindler & Diaferioll2008l - Chapter 17, this volume). 



3 Global abundances 

As a first step in our overview of the results of the ICM enrichment, we discuss a few general 
features of the enrichment pattern, which are related to the underlying model of chemical 
evolution. As a first example, we show in Fig.[3]maps of the fractional contribution of SN II 
to the global m etallicity obtained from SPH simulations of a galaxy cluster having a temper- 
ature c± 3 keV jTornatore et al.l2007ah : the simulations include a self-consistent description 
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Fig. 4 Iron mass ejection rate as a functi on of redshi ft for SNIa (top panel) and SNII (bottom panel), as 
predicted by the semi-analytical model by Cora (2006). The thick solid line is for the products released by 
the BCG, while the other curves are for galaxies lying within different cluster-centric distances, as computed 
at z = 0, after excluding the BCG contribution. 



of gas cooling, star formation and chemical evolution. The left and the right panels show the 
results based on assuming a Salpeter and a top-heavy AY IMFs, respectively. The clumpy 
aspect of these maps reveal that the products of SNII and SNIa are spatially segregated. 
Since SN II products are released over a relatively short time-scale, they are preferentially 
located around star forming regions. On the other hand, SNIa products are released over 
a longer time scale, which is determined by the life-times of the corresponding progeni- 
tors. As a consequence, a star particle has time to move away from the region, where it was 
formed, before SN la contribute to the enrichment. In this sense, the different spatial pat- 
tern of SN la and SN II p roducts also re flects the contribution that a diffuse population of 
intra-cluster stars (e.g., lArnaboldil l2004h for a review) can provide to the ICM enrichment. 
Also, using a top-heavier IMF turns into a larger number of SN II (see also Fig. 0, with a 
resulting increase of their relative contribution to the metal production. 

To further illustrate the different timing of the enrichment from SN la and SN II, we show 
in Fig.|4]the rate of the iron mass ejection from the two SN populations, as predicted by a 
SAM model of galaxy formation, coupled to a non-radiative hydrodynamical simulation 
dCora|2006). As expected, the contribution from SN II (lower panel) peaks at higher redshift 
than that of the SN la, with a quicker decline at low z, and closely traces the star formation 
rate (SFR). As for the SNIa, their rate of enrichment is generally given by the convolution 
of the SFR with the lifetime function. As such, it is generally too crude an approximation to 
compute the enrichment rate by assuming a fixed delay time with respect to the SFR. Also 
interesting to note from this figure is that everything takes places earlier in the BCG (solid 
curves), than in galaxies at larger cluster-centric distances. 



13 



A crucial issue when performing any such comparison between observations and sim- 
ulations concerns the definition of metallicity. From observational data, the metallicity is 
computed through a spectral fitting procedure, by measuring the equivalent width of an 
emission line associated to a transition between two heavily ionised states of a given ele- 
ment. In this way, one expects that the central cluster regions, which are characterised by a 
stronger emissivity, provides a dominant contribution to the global spectrum and, therefore, 
to the observed ICM metallicity. The simplest proxy to this spectroscopic measure of the 
ICM metallicity is, therefore, to use the emission-weighted definition of metallicity: 

_ ZiZiiniPgjAjTi.Zi) 

Y.i m iPg,iA\ T i) 

In the above equation, Z,, m,, p g j and 7] are the metallicity, mass, density and tempera- 
ture of the i-gas element, with the sum being performed over all the gas elements lying 
within the cluster extraction region. Furthermore, A (T,Z) is the metallicity dependent cool- 
ing function, which is defined so that n e n\\A(J \Z) is the radiated energy per unit time and 
unit volume fro m a gas element having elec tron number density « e and hydrogen number 
density «h (e.g. JSutherland & Dopitalll993l) . Since both simulated and observed metallic- 
ity radial profiles are characterised by significant negative gradients, we expect the "true" 
mass-weighted metallicity, 

I;Z;'M; 



(ID 



to be lower than the "observed" emission-weighted estimate. Indeed, mock spectral observa- 
tions of simulated clusters have shown that the above emis sion-weighted definition of metal- 
licity is generally quite close to the spectroscopic value dKapferer et al.ll2007al : iRasia et all 
2007), at least for iron. On the other hand, Rasia et al. have also shown that abundances of 
other elements may be significantly biased. This is especially true for those elements, like 
oxygen, whose abundance is measured from transitions taking place at energies at which it 
is difficult to precisely estimate the level of the continuum in a hot system, due to the limited 
spectral resolution offered by the CCD on-board Chandra and XMM-Ne wton . 

Based on a semi-analytical model of galaxy formation. iNagashima et all I J2005I) carried 
out a comparison between the global observed and predicted ICM enrichment for different 
elements. In their analysis, these authors considered both the case of a "standard" Salpeter 
IMF and the case in which a top-heavier IMF is used in starburst s triggered by galaxy m erg- 
ers. The comparison with observations, as reported in Fig. [5] led lNagashima et al] d2005h to 
conclude that the top-heavier IMF is required to reproduce the level of enrichment traced 
by different elements. All the abundances are predicted to be almost independent of tem- 
perature, at least for 7x > 1 keV (this temperature corresponds to the limiting halo circular 
velocity, above which gas and metals ejected by superwinds are recaptured by the halo). 
While this result agrees with the observational points shown in the figure, they are at vari- 
ance with respect to the s ignificant temperature dependence of the iron abundance found by 
iBaumgartner et al.| { |2b"05h . 

A different conclusion on the shape of the IMF is instead reached by other authors. 
For instance, Fabjan et al. (in preparation) analysed a set of hydrodynamical simulations 
of galaxy clusters, which include a chemical evolution model. In Fig. [6] we show the re- 
sult of this analysis for the comparison between the simulated temperature dependence of 
iro n (left panel) and s ulphur (righ t panel) and the results of XMM-Newton observations 
bv ldePlaaetalJ J20071) (see also I Werner et al.ll2008l, - Chapter 16, this volume). As for the 
iron abundance, there is a tendency for simulations with a Salpeter IMF to overproduce this 
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Fig. 5 The metal abundance s of the ICM, as predicted by a semi-analytical model of galaxy formation 
(from Nagashima et al. 2005). Each panel shows a different element: (a) [O/H], (b) [Fe/H], (c) [Mg/H] and 
(d) [Si/H]. Model predictions are shown by dots for a top-heavy IMF, and by crosses f or the Kennicutt 
IMF. Symbols with error bars correspond to observatio nal data, which have been take n from De Grand i et all 
2004), iPeterson etal] 120031) . iFukazawa eTal] <1998l) and iBaumgartner et alj )2005l) (seelNagashima et alj 



2005h for details). All abundances have been rescaled to the Solar abundances by Grevesse & Sauval 1 1998 



element, by about 30 per cent. Using a AY top-heavy IMF turns into a much worse dis- 
agreement, with an excess of iron in simulations by a factor 2-3. It is clear that this iron 
overabundance could be mitigated in case simulations include some mechanism for diffu- 
sion and/or mixing of metals, which decreases the abundance at small cluster-centric radii. 
Of course, this decrease should be such to reproduce the observed metallicity profiles (see 
below). As for the [S/Fe] relative abundancfl the difference between the two IMFs is much 
smaller, since the two elements are produced in similar proportions by different stellar pop- 
ulations, and, within the relatively large observational uncertainties, they agree with data. 

The comparison between the results of Figs.[5]and|6]is only one example of the different 
conclusions that different authors reach on the IMF shape from the study of the enrich- 



Here and in the following, we follow the standard bracket notation for the relative abundance of elements 
A and B: [A/B] = log /„ - log(Z A , Q /Z s , ). 
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Fig. 6 Left panel: the iron abundance within 0.2r5oo as a function of the cluster temperature (here rsoo is 
defined as the radius encompassing an overdensity of 500p c , with p c the cosmic critical density). The open 
squares are for simulated clusters analysed by Fabjan et al. (in preparation), which assume a Salpeter IMF, 
while the filled triangles are for two simulated clusters which assume the top-heavy AY IMF. Open and 
filled circle s with error bars are for cool- core and non cool- core clusters, respectively, observed with XMM- 
Newton by de Plaa et al. 1 200')}) (see also lWerner et al.l2 008 - Chapter 16, this volume). Right panel: the same 
as in the left panel, but for the abundance of sulph ur relative to iron. All abundances have been rescaled to 
the Solar abundances by Grevesse & Sauval 1 1998). 



ment of the ICM. Tracing back the origin of these differences is not an easy task. Besides 
the different descriptions of the relevant physical processes treated in the SAMs and in the 
chemo-dynamical simulations, a proper comparison would also require using exactly the 
same model of chemical evolution and the same sets of yields for different stellar popu- 
lations. We do recommend that all papers aimed at modelling the ICM enrichment should 
carefully include a complete description of the adopted chemical evolution model. 



4 Metallicity profiles 

Obse rvations based on ASCA (e.g. Fin oguenov et alj200ll) and Beppo-SAX (De Gran di et al.1 
|2004) satellites have revealed for the first time the presence of significant gradients in the 
profiles of ICM metal abundances. The much improved sensitivity of the XMM-Newton and 
Chandra satellites are now providing much more detailed information on the spatial distri- 
butions of different metals, thus opening the possibility of shedding light on the past history 
of star formation in cluster galaxies and on the gas-dynamical processes which determine 
the metal distribution. 

A detailed study of the metal distribution in the ICM clearly requires resorting to the 
detailed chemo-dynamical approach offered by hydrodynamical simulations. In principle, a 
number of processes, acting at different scales, ar e expected to play a role in determining th e 



spatial distribution of metals: galactic winds (e.g . Tomatore et al. 2004; Romeo et alj200d) 
transport by buoyant bubbles ([Roediger et al. 12007 ; Siiacki et alj 2007 ), ram-pressure strip - 
ping (e.g jKapferer et al. 2007b), diffusion by stochasti c gas motions dRebusco et al . 2005) 



sinking of highly enriched low-entropy gas dCoral r2006). etc. It is clear that modelling prop- 
erly all such processes represents a non-trivial task for numerical simulations, which are 
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Fig. 7 Left panel: contribution from ram-pre ssure stripping and f rom galactic winds to the resulting iron 
metallicity profile of a simul ated cluster (f rom Kapferer et al. 2007b). Data points are from the Beppo-SAX 
observations analysed by De Grandi et al. (2004). Right panel: the effect of AGN feedback on t he profiles of 
mass- weighted gas metallicity and stellar metallicity with simulations of a galaxy cluster (from Siiac ki et al.l 
2007). Continuous lines show the metallicity profiles without AGN heating, while dashed lines are for the 
runs with AGN feedback. The arrows illustrate how the ICM metallicity is affected by the AGN feedback. 
The vertical dotted lines mark the resolution limit of the simulations and the virial radius of the simulated 
cluster. 



required both to cover a wide dynamical ra nge and to provide a reliabl e description of a 
number of complex astrophysical processes dSchindler & Diaferioll2008l - Chapter 17, this 
volume). 

As an example, in the left panel of Fig.[7]we show the result of the analysis by Kapferer et al. 
J2007bl) . which was aimed at studying the relative role played by galactic winds and ram- 
pressure stripping in determining the metallicity profiles. In these simulations, based on an 
Eulerian grid code, these authors use a semi-analytical model to follow the formation and 
evolution of galaxies, and described the evolution of a global metallicity, instead of follow- 
ing different metal species. It is interesting to see how different mechani sms contribute in the 
simul ations to establish the shape of the metallicity profile (see also lSchindler & Diaferiol 
120081 . - Chapter 17, this volume). In particular, ram-pressure stripping turns out to be more 
important in the central regions, where the pressure of the ICM is higher. On the other hand, 
galactic winds are more effective at larger radii, where they can travel for larger distances, 
thanks to the lower ICM pressure. 

As a further example, we show in the right panel of Fig.|7]the effect of AG N feedback on 
the di stribution of metals in both gas and stars from the cluster simulations bv lSiiacki et alj 
d2007l) . These authors included in their SPH simulations a model for black hole growth and 
for the resulting energy feedback. Since these simulations do not include a detailed model 
of chemical evolution, the resulting metallicity cannot be directly compared to observations. 
Still, this figure nicely illustrates how the effect of AGN feedback is that of changing the 
share of metals between the stars and the gas. Due to the quenching of low-redshift star 
formation, a smaller amount of metals are locked back in stars, with a subsequent reduction 
of their metallicity. Furthermore, the more efficient metal transport associated to the AGN 
feedback causes a significant increase of the ICM metallicity in the outer cluster regions, 
thereby making the abundance profile shallower. 
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Since iron is the first element for which reliable profiles have been measured from X-ray 
observations, most of the comparisons performed so far with simulations have concentrated 
on its spatial distribution. As an exam ple, we show in the left panel of Fig.[8]the comparison 
between the iron profiles meas ured by De Grandi et al.l fc004l) and results from the chemo- 
dynamical SPH simulations by iRomeo et alJ d200fjh . These authors performed simulations 
for two clusters having sizes comparable to those of the Virgo and of the Coma cluster, 
respectively. After exploring the effect of varying the feedback efficiency, they concluded 
that the resulting metallicity pro files are always st eeper than the observed ones. A similar 
conclusion was also reached by IValdarninU ( l2003h . based on SPH simulations of a larger 
ensemble of clusters, which used a rather inefficient form of thermal stellar feedback. A 
closer agreement between observed and simulated profiles of th e iron abundance is found 
by Fabjan et al. (in preparation, see also iTornatore et al]|2007ah . based on GADGET-2 SPH 
simulations, which include the effect of galactic winds, pow ered by SN explosions, accord- 
ing to the scheme presented bv lSpringel & Hernquistl J2003n . The results of this comparison 
are shown in the right panel of Fig. [8] The solid black lines show the results for simulated 
clusters of mass in the range ~ (1 — 2) x 10 15 h~ 1 M ( j ) , based on a Salpeter IMF, while the 
dashed curve is for one cluster simulated with a top-heavy AY IMF. The observational data 
points are for a set of nearby relaxed clusters observed with Chandra and for two clusters 
with T c± 3 keV observed with XMM-Newton. This comparison shows a reasonable agree- 
ment between simulated and observed profiles, with the XMM-Newton profiles staying on 
the low side of the Chandra profiles, also with a shallower slope in the central regions. While 
this conclusion holds when using a Salpeter IMF, a top-heavy IMF is confirmed to provide 
too high a level of enrichment. 

The different degree of success of the simulations shown in the two panels of Fig. [8] in 
reproducing observations can have different origins. The most important reason is probably 
due to the different implementations of feedback which, as shown in Fig. [7] plays a deci- 
sive role in transporting metals away from star-forming regions. Another possible source 
of difference may lie in the prescription with which metals are distributed around star par- 
ticles. Once again, the details of the numerical implementations of different effects need 
to be clearly specified in order to make it possible to perform a meaningful comparison 
among different simulations. Finally, it is also interesting to note that clusters observed with 
Chandra seem to have steeper central profiles than those observed with XMM-Newton and 
Beppo-SAX. There is no doubt that increasing the number of clusters with reliably deter- 
mined metallicity profiles will help sorting out the origin of such differences. 



5 Evolution of the ICM metallicity 

After the pioneering paper based on ASCA data by Mushot zkv & Loewensteinl l ll997h . the 
increased sensitivity of the XMM-Newton and Chandra observatories recently opened the 
possibility of tracing the evolution of the global iron c ontent of the ICM out to the largest 
redsh ifts, z — 1.4, where clusters have been secured (e.g. Balestra et al. 2007; Maughan et al. 
120071) . These observations have shown that the ICM iron metallicity decreases by about 50 
per cent from the nearby Universe to z — 1 • This increase of metallicity at relatively recent 
cosmic epochs can be interpreted as being due either to a fresh production of metals or to a 
suitable redistribution of earlier-produced metals. Since SN la, which provide a large contri- 
bution to the total iron budget, can release metals over fairly long time scales, the increase 
of the iron abundance at z < 1 is not necessarily in contradiction with the l ack of recent 
star formation in the bulk of cluster galaxies. On the one hand, recent analyses tettorill2005l ; 
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Fig. 8 Left panel: profiles of iron metallicity at z = for different runs of two simulated clusters, obtained by 
changing the IMF (AY: Arimoto and Yoshii 1987; Sal: Salpeter 19 55), the efficiency of the energy feedback 
and by including the effect of thermal conduction (COND; from Romeo et al. 2005). Observational data 
points are taken from De Grandi et al. (2004). Right panel: Comparison between observed and simulated 
profiles of iron abundance (from Fabjan et al., in preparation). The continuous black curves are for four 
simulated clusters, having virial masses in the range (1 — 2) X 1O 15 /i~'M0, assuming a Salpeter IMF. The 
black dashed line is for one of these clusters, simulated b y assuming the AY top-heavy IMF. Filled circles 
with error bars shows observational results from Chandra I Vikhlinin et al. 2005) for clusters with temperature 
above 3 keV, whi le the open square s and triangles are f or the profiles of two ~ 3 keV clusters observed with 
XMM-Newton i de Plaa et al. 2006; Werner et al. 2006). The iron abundance has been rescaled to the Solar 
abundances by Grevesse & SauvaJ dl998l) . 



iLoewensteinl l2006h have shown th at the observed metal licity evolution is consistent with 
the observed SN la rate in clusters -Yam et al. I l2002h . On the other hand, a number of 
proce sses, such as stripping of metal enriched gas from merging galaxies (e.g. JCalura et alj 
l2007h or sin king of highly enriched low-entropy gas towards the cluster central regions 
i Con]|2006l) may also significantly contribute to the ICM metallicity evolution, through a 
redistribution of enriched gas. Indeed, the amount of accreted gas since z — 1 is at least 
comparable to that present within the cluster at that redshift. Therefore, the evolution of the 
ICM metallicity is expected to be quite s ensitive to the e nrichment level of the recently ac- 
creted gas. On the other hand, as noted bv lRenzinil ( ll997l) , if ram-pressure stripping plays a 
dominant role in determining the ICM enrichment, one should then expect to observe more 
massive clusters, whose ICM has a higher pressure, to be more metal rich than poorer sys- 
tems. In fact, observations do not show thi s trend and, if any, sugges t a decrease of the iron 
abundance with the ICM temperature (e.g.. lBaumgartner et alj|2005h . 

ICoral ( f2006h used a hybrid approach based on coupling a SAM model of galaxy forma- 
tion to a non-radiative hydrodynamical simulation to study the gas dynamical effects which 
lead to the re-distribution of metals from the star forming regions. The results of her analysis 
are summarised in the metallicity maps shown in Fig. [9] Metals are produced inside small 
halos which at high redshift define the proto-cluster region. Since the gas in these halos has 
undergone no significant shocks by diffuse accretion, its entropy is generally quite low. This 
low-entropy highly enriched gas is therefore capable to sink towards the central regions of 
the forming cluster. As a consequence, this leads to an increase of the metallicity in the re- 
gions which mostly contribute to the emissivity, thereby causing a positive evolution of the 
observed metallicity. 
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Fig. 9 Evolution of the iron abundance within 700 h kpc from the centre of a simulated cluster at 
z = 1,0.5,0.3,0.2,0.1,0, from upper left to bottom right panel. The maps show the projection of the mass- 
weighted ir on abundance by number relative to hydrogen, Fe/H, with respect to the solar value (adapted from 

ICo^l2006l) . 



Nagashima et al. U2005h used a SAM approach, not coupled to hydrodynamical simula- 
tions, to predict the evolution of the ICM metallicity. In this case, the evolution is essentially 
driven by the low-redshift star formation, since any redistribution of enriched gas is not 
followed by the model. The result of their analysis (left panel of Fig. 1 10b shows that a very 
mild evolution is predicted. A prediction of the evolution of the iron abundance from hy- 
drodynamical simulations, which self-consistently follow gas cooling, star formation and 
chemical evolution, has been more recently performed by Fabjan et al. (in preparation). The 
results of their analysis are sh own in the right panel of Fig. [TO] where they are also compared 
to the observational results by Bale stra et al.l d2007[) . Also in this case, the results based on a 
Salpeter IMF are in reasonable agreement with observations, while a top-heavier IMF pro- 
duces too high metallicities at all redshifts. Quite interestingly, at low redshift the run with 
top-heavy IMF shows an inversion of the evolution. This is due to the exceedingly high star 
formation triggered by heavily enriched gas, which locks back in stars a substantial fraction 
of the gas metal content in the cluster central regions. This highlights the crucial role played 
by star formation in both enriching the ICM and in locking back metals in the stellar phase. 
This further illustrates how the resulting metallicity evolution must be seen as the result of 
a delicate interplay of a number of processes, which all need to be properly described in a 
self-consistent chemo-dynamical model of the ICM. 
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Fig. 10 The redshift evolution of the iron abundance. Left panel: predictions by different flavours of the SAM 
by Nagashima et al. (2005), based on a top-heavy IMF, compared with observational data from different 
sources, for clusters within two temperature ranges. Error bars on model predictions are the lo~ scatter over 
the ensemble of cluster formation histories. Right panel: comparison between observational data and full 
hydrodynamical si mulations by Fab j an et al. (2007, in preparation). Open circles come from the analysis 
of Chandra data by Balestra et al. (2007). The shaded area represents the r.m.s. scatter within each redshift 
interval, while the error bars are the 1 o~ confidence level. The filled squares are for a set of four simulated 
clusters, using a Salpeter IMF, with temperature kT > 5 keV, with error bars indicating the corresponding 1 o~ 
scatter. The filled triangles refer to only one cluster, out of four, simulated with the top-h eavy IMF. In both 
panels , the values of the iron abundance have been rescaled to the Solar abundances by Grevesse & Sauval 



6 Properties of the galaxy population 

The physical processes which determine the thermodynamical and chemical properties of 
the ICM are inextricably linked to those determining the pattern of star formation in cluster 
galaxies. For this reason, a successful description of galaxy clusters requires a multiwave- 
length approach, which is able at the same time to account for both the X-ray properties 
of the diffuse hot gas and of the optical/near-IR properties of the galaxy populations. In 
this sense, the reliability of a numerical model of the ICM chemical enrichment is strictly 
linked to its ability to reproduce the basic properties of the galaxy population, such as lumi- 
nosity function and colour-magnitude diagram. While the study of the global properties of 
the galaxy popula tions has been conducted for several years wi th semi-analytical models of 
galaxy formation jDe Lucia et al.l20 04; Nagashim a et al.l2005h . it is only quite recently that 
similar studies have been performed with cosmological hydrodynamical simulations, thanks 
to the recent advances in supercomputing capabilities and code efficiency. 

As we have already mentioned, hydrodynamical simulations treat the process of star for- 
mation through the conversion of cold and dense gas particles into collisionless star particles. 
Each star particle is then characterised by a formation epoch and a metallicity. Therefore, it 
can be treated as a SSP, for which luminosities in different bands can be computed by resort- 
ing to suitable stellar population synthesis models. Since the colours of a SSP are quite sen- 
sitive to metallicity, it is necessary for simulations to include a detailed chemical evolution 
model in order to reliably compare predictions on optical/near-IR properties of galaxies with 
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Fig. 11 Left panel: Colour-magnitude relations for the simulations of elliptical galaxies I Kawata & Gibson 
2003b), based on a Salpeter IMF, compared to the observational results for the galaxies in Coma ( Bower et al. 
1992) (crosses) and for NGC 4472 (open star). The filled square and circle are the results for the reference 
runs, whereas the open symbols are when the contribution from stars younger than 2 or 8 Gyr are ignored. 
Right panel: the colour-magnitude relation for simulated clusters (symbols). Upper an d lower panels re fer to 
simulations performed with a Salpeter IMF and with the AY IMF, respectively (from Saro et tnriooa) . The 
big dots are for the BCGs of the simulated clusters. The straight lines indicate the best-fitting relation and the 
scatter observed for the Coma galaxies. 



observational data. This approach has been pursued by Kawata & Gibson! (2003b), who per- 
formed simulations of elliptical galaxies and bv lRomeo et al.l ( 120051 ) and ISaroet alJ (12006). 
who performed simulations of galaxy clusters. The results of these analyses are summarised 
in Fig. [TT] In the left panels we show the position on the c olour-magnitude diagram of the 
elliptical galaxies simulated bv iKawata & Gibson! J2003bh . This plot clearly indicates that 
elliptical galaxies are generally much bluer than observed, as a result of an excess of recent 
star formation, which is not properly balanced by an efficient feedback mechanism. Indeed, 
it is onl y after excluding the contribution from recently formed stars that lKawata & Gibson! 
(2003b) were able to obtain ellipticals with more realistic colours. 

In the right panels we show the colour-magni tude relation (C MR) of the galaxy popu- 
lation of an extended set of clusters simulated bv lSaro et al.l ([2006) based both on a Salpeter 
IMF and on a top-heavy IMF. Quite apparently, the bulk of the galaxy population traces 
a CMR which is consistent with the observed one, at least for a Salpeter IMF. As for the 
top-heavy IMF, it predicts too red colours as a consequence of a too high metal content of 
galaxies. Also interesting to note, brighter redder galaxies are on average more metal rich, 
there fore demonstrating that the CMR corresponds in fact to a metallicity sequence (see 
also iRomeo etal 1 12003) . However, a clear failure of these simulations is that the BCGs are 
always much bluer than observed. This result agrees with the exceedingly blue colour of 
simulated ellipticals shown in the left panels. This is due to an excess of recent star forma- 
tion, which, in B CGs of ~ IQ^Mp, s imulated clusters, can be as high as 500-1000 M yr~' . 
As discussed by Borganietal ] |2008l - Chapter 13, this volume, the presence of overcooling 
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in cluster simulations manifests itself with the presence of steep negative temperature gra- 
dients in the central regions. In this sense, wrong BCG colours and temperature profiles in 
core regions are two aspects of the same problem, i.e. the difficulty of balancing the cooling 
runaway within simulated clusters. 

Since all these simulations include recipes for stellar feedback, although with different 
implementations and efficiencies, the solution to the overcooling in clusters should be found 
in some non-stellar feedback, such as that provided by AGN. Recent studies of the galaxy 
properties, based on semi-analytical models, have shown that an energy feedback of the 
sort of that expected from AGN, is able to suppress the bright end of the galaxy luminosity 
function, which is populated by the b ig ellipticals residing at the centre of clusters (e.g., 
iBower et alj|2006l : ICroton et alj|200d) . There is no doubt that providing robust predictions 
on the properties of the galaxy population will be one of the most important challenges for 
simulations of the next generation to be able of convincingly describing the thermo- and 
chemo-dynamical properties of the hot gas. 



7 Summary 

The study of the enrichment pattern of the ICM and its evolution provides a unique means 
to trace the past history of star formation and the gas-dynamical processes, which deter- 
mine the evolution of the cosmic baryons. In this overview, we discussed the concept of 
model of chemical evolution, which represents the pillar of any chemo-dynamical model, 
and presented the results obtained so far in the literature, based on different approaches. Re- 
stricting the discussion to those methods which follow the chemical enrichment during the 
cosmological assembly of galaxies and clusters of galaxies, they can be summarised in the 
following categories. 

(1) Se mi-analytical models (SAMs) of galaxy formation dPe Lucia et al.l2004l ; lNagashima et all 
l2005h : the production of metals is traced by the history of galaxy formation within the evolv- 
ing DM halos and no explicit gas-dynamical description of the ICM is included. These ap- 
proaches are computationally very cheap and provide a flexible tool to explore the space 

of parameters determining galaxy formation and chemical evolution. A limitation of these 
approaches is that, while they provide predictions on the global metal content of the ICM, 
the absence of any explicit gas dynamical description causes the lack of useful information 
on the spatial distribution of metals. 

(2) SAMs coupled to hydrodynamical simulations. In this case, the galaxy formation is fol- 
lowed as in the previous approach, but the coupling with a hydrodynamical simulation al- 
lows one to trace the fate of the enriched gas and, therefore, to study the spatial distribution 
of metals. Different authors have applied different implementati ons of this h ybrid approach, 
by focusing either on the role of the chemical evolution model l|Cora|[2006) or on the effect 
of specific gas-dynamic al processes, such as r am-pressure stripping and galactic winds, in 
enriching the ICM (e.g., Ka pferer et al.ll2007bh . The advantage of this approach is clearly 
that gas-dynamical processes are now included at some level. However, the galaxy forma- 
tion process is not followed in a self-consistent way from the cooling of the gas during the 
cosmic evolution. 

(3) Full hydrodynamical simulat ions, which self- c onsistently follow gas cooling, star for- 
mation and chemical evolution JValdarninill2003l ; iKawata & Gibsonl l2003bl : iRomeo et alj 
l2005l ; lTornatore et al.l l2007ah . Rather than being described by an external recipe, galaxy for- 
mation is now the result of the cooling and of the conversion of cold dense gas into stars, 
as implemented in the simulation code. The major advantage of this approach is that the 
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enrichment process now comes as the result of the star formation predicted by the hydro- 
dynamical simulation, without resorting to any external model. The bottleneck, in this case, 
is represented by the computational cost, which becomes prohibitively high if one wants to 
resolve the whole dynamic range covered by the processes of metal production and distribu- 
tion, feedback energy release, etc. For this reason, these simulation codes also need to resort 
to sub-resolution models, which provide an effective description of a number of physical 
processes. However, there is no doubt that the ever improving supercomputing technology 
and code efficiency make, in perspective, this approach the way to study the past history of 
cosmic baryons. 

As already emphasised, developing a reliable model of the cosmic history of metal pro- 
duction is a non-trivial task. This is due to the sensitive dependence of the model predictions 
on both the parameters defining the chemical evolution model (i.e., IMF, stellar life-times, 
stellar yields, etc.), and on the processes, such as ejecta from SN explosions and AGN, 
stripping, etc., which at the same time regulate star formation and transport metals outside 
galaxies. In the light of these difficulties, it is quite reassuring that the results of the most 
recent hydrodynamical simulations are in reasonable agreement with the most recent obser- 
vational data from Chandra and XMM-Newton. 

In spite of this success, cluster simulations of the present generation still suffer from 
important shortcomings. The most important of them is probably represented by the diffi- 
culty in regulating overcooling at the cluster centre, so as to reproduce at the same time the 
observed "cool cores" and the passively evolving massive ellipticals. Besides improving the 
description of the relevant astrophysical processes in simulation codes, another important 
issue concerns understanding the possible observational biases which complicate any direc t 
comparison between data and model predictions dKapferer et afl|2007al : iRasia et alj[2007h . 
In this respect, simulations provide a potentially ideal tool to understand these biases. Mock 
X-ray observations of simulated clusters, which include the effect of instrumental response, 
can be analysed exactly in the same way as real observational data. The resulting "observed" 
metallicity can then be compared with the true one to calibrate out possible systematics. 
There is no doubt that observations and simulations of the chemical enrichment of the ICM 
should go hand in hand in order to fully exploit the wealth of information provided by X-ray 
telescopes of the present and of the next generation. 
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